International Journal of Electrical and Computer Engineering (IJECE) 
Vol. 12, No. 2, April 2022, pp. 1795~1806 
ISSN: 2088-8708, DOI: 10.1159 1/ijece.v1212.pp1795-1806 O 1795 


A hybrid approach for categorizing images based on complex 
networks and neural networks 


Ali Ebrahimi!, Kamal Mirzaie!, Ali Mohamad Latif? 
'Department of Computer Engineering, Maybod Branch, Islamic Azad University, Maybod, Iran 
Department of Computer Engineering, Yazd University, Yazd, Iran 


Article Info ABSTRACT 

Article history: There are several methods for categorizing images, the most of which are 
; statistical, geometric, model-based and structural methods. In this paper, a 

Received Mar 13, 2021 new method for describing images based on complex network models is 

Revised May 31, 2021 presented. Each image contains a number of key points that can be identified 

Accepted Jul 1, 2021 through standard edge detection algorithms. To understand each image 


better, we can use these points to create a graph of the image. In order to 
facilitate the use of graphs, generated graphs are created in the form of a 
Keywords: complex network of small-worlds. Complex grid features such as topological 
and dynamic features can be used to display image-related features. After 
generating this information, it normalizes them and uses them as suitable 


Complex network 


Key points ee features for categorizing images. For this purpose, the generated information 
Network similarity is given to the neural network. Based on these features and the use of neural 
Neural network networks, comparisons between new images are performed. The results of 
Small-world network the article show that this method has a good performance in identifying 
similarities and finally categorizing them. 

This is an open access article under the CC BY-SA license. 
[OK] 

Corresponding Author: 


Kamal Mirzaie 

Department of Computer Engineering, Maybod Branch, Islamic Azad University 
Maybod, Iran 

Email: k.mirzaie @maybodiau.ac.ir 


1. INTRODUCTION 

With the production of the first digital image in 1957 by Russell Kirsch a scientist at the National 
Research Institute for Standards and Technology (NIST) the world has virtually entered a new era of digital 
imaging. Rapid development of digital imaging intensified after the advent of microprocessors in the early 
1970s. Millions of digital imaging machines have revolutionized conventional photography over the past 
decade and into the late twentieth century, producing billions of digital images. Due to the huge volume of 
digital images on Internet and their diversity, large and small digital libraries have been produced around the 
world, which has led to manage image databases and use tools for image retrieval and classification [1]. 

The term "content-based image retrieval" dates back to 1992 when Kato used a database based on 
existing colors and shapes to describe experiments involving automatic image retrieval. Since then, this term 
has been used to describe the process of retrieving images from a large set based on coordinated visual 
features. The used techniques, tools, and algorithms apply methods such as statistics, pattern recognition, 
signal processing, and computer vision. The first commercial content-based image retrieve (CBIR) system 
developed by immersed boundary method (IBM) was called query by image content (QBIC) [2]. 

Textual information about images can be easily searched using existing technologies, but it is a 
human-dependent method to manually describe the information of each image in the database. Therefore, this 
method is practically not applicable for very large databases or images that are generated automatically. The 
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original CBIR system was designed to search databases based on color, texture, and shape properties. After 
the development of these systems, the need for a user-friendly interface became apparent. Therefore, CBIR 
includes a plan that tries to meet the needs of the user and perform his searches. Conventional methods for 
image retrieval include: semantic evaluation, feedback method (two-way communication between machine 
and human), machine learning, query methods, content comparison using image distance measurement, 
texture, correlation matrix, and wavelet transform are all used [2]. 

So far, different methods with advantages and disadvantages have been introduced to identify image 
patterns. In recent years, the recognition of texture and the classification of images by color and form have 
been in the center of attention. For example, Ramesh introduced a histogram-based system that uses 
histogram characteristics to categorize images [3]. Moreover, Shirazi et al. [1] compared the various methods 
used by the vector machine for categorizing the images. Shakoor et al. [2] used a local entropy method to 
identify the image features and used them for categorization. 

As mentioned, one of the most important features of any image is its key points. In other words, key 
points are one of the most important indicators of any image, and so far, various algorithms have been 
designed to identify them, some of which are more important or are used more. Edges can be used as one of 
the key points in identifying and classifying images. The purpose of identifying edges is to mark points in an 
image where the intensity of light changes dramatically. Rapid changes in image properties usually represent 
important events and changes in environmental characteristics [3]. 

The basic theory in most edge detection methods is the calculation of a local derivative operator. 
The size of the first derivative of the image at each point is equal to the size of the gradient. The second 
derivative is also obtained using the Laplace approximation. If we consider an edge as a change in brightness 
seen over several pixels, edge detection algorithms generally calculate a derivative of this change in 
brightness. After calculating the derivative, the next step is to apply a threshold to discover the points that are 
part of an edge. The lower the threshold, the more lines are detected and the more sensitive the results are to 
noise and unrelated image features. On the other hand, a high threshold may miss weak lines or parts of lines. 
In most works done to retrieve images or determine the similarity between two images, edge recognition has 
been used to identify key points. For this purpose, edge finding algorithms such as Harris, Shi, Sift, and Surf. 
have been used. As the name implies, these algorithms are used to identify edges or corners of an image that 
are more important than other points. Each of these algorithms has been used in specific cases, but in general, 
each has been created to complete the shortcomings of other existing algorithms and to create a more 
accurate algorithm for edge and key point detection, although a complete and accurate edge algorithm for all 
images has not been designed yet [4]. 

For example, the scale invariant feature transform (SIFT) algorithm is used in tasks such as 
matching different faces of a body or scene, and identifying objects. The characteristics obtained are 
somewhat constant regardless of the image scale, rotation, change of view, and lighting variations [4]. The 
speed-up robust feature (SURF) algorithm is also a feature and descriptor detection method, and its standard 
version is much faster than SIFT. The algorithm with an accurate and decisive approximation is used to 
detect points, which can be calculated with the help of the integral. SURF is used to detect objects, 
individuals, or faces to rebuild 3D scenes for tracking objects [5]. Another algorithm is the Shi-Tomasi 
method, in which the minimum of the values is computed. Based on the matrix theory, these values are used 
to identify the edges and key points [6]. One of the other algorithms is the Harris method, which is 
introduced by Harris and Stephens and improved the moravec corner detector method [7]. 

Tomasi and Kanade [6], founded the theory of graphs in 1736 by solving the problem of the 
Konigsberg Stairs. But James Joseph Sylvester was the first to call these mathematical models graphs in 
1878. Graphs play an important role in simplifying images and their understanding by computer systems. As 
we know, each graph consists of a number of nodes that are connected to each other by edges. Graphs have 
different models (directed graph, undirected graph, mixed graph, weighted graph, oriented graph, and cycle 
graph). 

Graph similarity determination is particularly challenging as graphs fail to be naturally represented 
or converted into vector spaces. Likewise, graphs intended for pattern recognition have many edges and 
nodes, rendering them are busy. Graph similarity can be measured in a variety of ways. Main methods are 
based on graph spacing in graph edit distance (GED) [8], [9], the graph core [10]-[12], and graph matching 
and topological description methods [12]. Sanfeliu and Fu [8] pioneered the notion of GED, which is being 
widely used. Graphs can be typically converted to string sequences, with their direction utilized to compare 
these strings [13]. Using the Levenshtein distance, Hancock calculated string similarity [14]. Tang converted 
the graphs to the string, and based on the theory of the catch graph, and using the dynamic time deviation of 
dynamic time warping (DTW), calculated the similarity of the fibers [15]. This problem was formulated by 
Kelly and Hancock [13] to increase the likelihood of series-mediated binary graph alignment. The Gao 
algorithm distance (GAD) provided a histogram based on GED calculation for two graphs [16]. 
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Graph similarity can also be measured by calculating the similarity between a couples of frequently- 
used pattern-based graphs. There are simple as well as complex patterns. Core graph function development 
practices have been thoroughly examined in the best possible way. These methods can be broadly categorized 
based on sub-classes [8], shortest paths [5], random steps [15] and sub-graphs [12]. 

In stochastic nuclei, graph similarity is measured by counting basic steps between each pair. This 
kind of core is usually very complex and shaky. By replacing the steps with paths, Bargwardt et al. 
introduced the shortest nucleus of the path on the graphs [11]. The main goal is to map the graphs into vector 
spaces in which important dimensions are concerned. Furthermore, various image processing methods can be 
investigated by reviewing the manifold method if a given graph can be embedded in manifold template 
space. Nevertheless, a number of counterparts are applicable in creating low-dimensional spaces for a graph 
set. This can be accomplished by representing graphs using pattern vectors. This problem can be solved using 
topological descriptors [16], [17]. These methods first map out each graph and then use the spacing and 
metric on the vector to study the graphs. Many actions have been taken to illustrate the topology of the graph. 
For example, research [18], [19] mapped the structure of a graph into a fixed-length design through the 
spectral graph theory. They exploited vectors derived from graphs to manage the calculation of the distance 
between two graphs and evaluate their structures according to principal component analysis (PCA) and multi- 
dimensional scaling (MDS) techniques. 

The Gao graph with the histogram schematics [20] is similar to the Luo method [18]. In this method, 
the cost function or node similarity level in similar vectors in the two graphs cannot be formulated. The 
determination of both cases may be tricky, rendering the method disadvantageous mainly because it fails to 
examine graph structures in accordance with histogram plotters. The present article used the method of 
displaying complex networks. Complex networks are characterized as taking advantage of the combination of 
graph theory, probability theory. Recent efforts have led to the use of complex networks in various 
applications and fields. For example, Khani and Mirzaei used the structure of the small network for 
simulating cancer as a multi-functional approach to biological systems [20]. 

Liu et al. [21] provided a time series clustering given the high dimensions of multivariate time series 
and most previous work focusing on univariate time series clustering, a new method based on principal 
component analysis to achieve clustering. Multivariate time series are presented with more speed and 
accuracy. They are inspired by the traditional k-Means clustering method and have used the detailed 
algorithm of the proposed Mc2PCA method to analyze and cluster multivariate time series. 

Lü et al. [22] provided a new technology to research the problem of the synchronization 
transmission of spiral wave and turbulence between uncertain time-delay neuronal networks with different 
node numbers and topologies. Recently novel ideas about advancing evolutionary dynamics towards new 
phenomena including many new topics, even the dynamics of equivalent social networks have been 
proposed. In fact, they include more advanced complex networks and incorporate them with the coupled map 
lattices (CMLs), which are usually used for spatiotemporal complex systems simulation and analysis [23], 
[24]. 

This study graphically represented the structure of a complex network-based grid. The graph 
structure is initially mapped as a complex network, followed by the description of the graph structure by 
deriving the topological and dynamic complexity of the complex network. These features were extracted 
using neural networks and will be used to categorize images. 

The study of complex networks is a young and active area of scientific research (since 2000) 
inspired largely by empirical findings of real-world networks such as computer networks, biological 
networks, technological networks, brain networks, climate networks and social networks. Euler's solution of 
the Seven Bridges of Konigsberg problem is considered to be the first true proof in the theory of networks. 

Complex networks can be categorized into three models: i) random networks, ii) small-world 
networks, and iii) scale-free networks. In the random model, which is the simplest model, the edges are 
added randomly. A small-world network is a type of mathematical graph in which most nodes are not 
neighbors of one another, but the neighbors of any given node are likely to be neighbors of each other and 
most nodes can be reached from every other node by a small number of hops or steps. Specifically, a small- 
world network is defined to be a network where the typical distance L between two randomly chosen nodes 
(the number of steps required) grows proportionally to the logarithm of the number of nodes N in the 
network. 

A scale-free network is a type of random network in which some nodes have been removed. In this 
network, a small number of nodes are very important and by removing them, the network is broken. These 
nodes are called poles. Social networks are scale-free networks. 

The complex network edges as a graph can be binary, weight-oriented, and directional. In this 
article, only an undirected model was checked. The non-directional complex network model was usually 
represented by the W matrix, where w (i, j) is the weights that connect node i to node j. 
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The topology and complexity of the network can be detected using the measurements taken from the 
graph theory and also for complex network research. Every complicated network has a specific topological 
design that determines how to connect. Accordingly, complex networks can be transformed depending on 
what measures are used, facilitating the derivation and classification of suitable topological designs . 

Grade distribution is an important attribute of the vertex in the graph. Given vertex degrees, many 
measures can be formulated, with the two most straightforward measurement techniques as in (1) and (2): 


Kmax = max; K; (1) 
n 
K, 1 >: K, 
avg 77 i 
ni (2) 


where K; indicates the degree of node i. Structural and dynamic networks can be principally characterized by 
connectivity degree distribution, including: i) entropy, ii) energy, and iii) average connection degree. 
The entropy of connectivity degree distribution can be determined as shown in (3): 


Kmax 


Kalk’) =- Ù, p(k,k') log (rk, k’)) 
k,kr=1 (3) 


The normal distribution-associated energy can be determined as shown in (4): 


Kmax 
E= > pak’) 
k,kr=1 (4) 


One way to detect the loop is through the clustering coefficient, which is beneficial for network analysis like 
degree. Two different clustering coefficients are typically used, the first of which, according to the definition 
of nondirectional networks, is as in (5): 


3N; 
E = — 
N; (5) 
where Ny, represents the number of triangles on the grid and N; represents the triangles. 

A triangle has three vertices, with an edge between each pair, a triangular connected set of three 
vertices, whose two vertices are adjacent to a third one (i.e., central vertex). Clustering factors can be 
obtained for a grid by calculating the mean clustering factor for each vertex. The clustering factor of vertex i 
can be measured as in (6): 


es Na (i) 
‘Ns (i) (6) 


where N,(i) denotes the number of triangles with vertices i and N3(i) denotes the number of triangles, i.e., 
the central vertex. The clustering coefficient of a network can be defined as in (7): 
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The distance is an important size, depending on the overall network structure. We can measure the network 
as the mean of the shortest distance by calculating the mean distance of the shortest line (path) for each pair 
of vertices, such as (8): 
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(7) 


(8) 


dij represents the shortest line distance from vertex i and vertex j. 


Int J Elec & Comp Eng, Vol. 12, No. 2, April 2022: 1795-1806 


Int J Elec & Comp Eng ISSN: 2088-8708 O 1799 


2. DISPLAY OF COMPLEX NETWORK PROPERTIES 

This section initially models a graph by a complex network and, consequently, extracts the dynamic 
and topological properties of the above network to detect the graph structure, which is named “Display of 
complex network properties of a complex network properties extraction (CNCRG) graph”. Thus, the distance 
between CNCRGs can be taken as a similarity index between the two graphs. Regarding the relation between 
graph edges constituting the basic elements of the graph structure, CNCRGs would help identify the 
differences between the two graphs. This can be associated with GED, which takes into account differences 
in the binding of the graphs through the order of editing operations, including adding, deleting, and replacing 
the edges. The CNCRG is completely independent of the cost function. 


3. PROPOSED APPROACH BASED ON COMPLEX NETWORKS AND NEURAL NETWORKS 

In this paper, the proposed method was a combination of complex networks and neural networks. 
Using the complex network, the properties were extracted and then, the neural network was used to classify. 
Assume graph G=(V, E). Here, E and V indicate the set of edges and vertices, respectively. |E| and V| denote 
the number of edges and nodes of the graph G, respectively. The network model can be defined as 
GN=(Vy, En) for the G graph. Each edge of the graph G is a vertex GN. The weighting vector for each vertex 
can be measured by evaluating the relationship between and geometric position of the edges, as shown in 
Figure 1, the weighting vector for vertex i can be measured as seen in (9): 


e; = (l, di, dii dzi, Xi Yi) (9) 


where l; represents edge length, d; represents the distance between the edge center and graph center (Point o), 
dz; and d4; the intervals from the beginning of the edge to the center of the graph (point o) and from edge end 
to graph center (point o) and x; and y; are the center point coordination of the edges. 


Figure 1. Schematic representation of a graph based on a complex network 


The network is formed by a set of nondirectional edges E that connects each pair of vertices. The 
value of the edge of vertex i to vertex j can be calculated using the Euclidean distance between e; and ej: 


Wi; = d(e;, e;) (10) 
Therefore, the network can be mapped as |E|x|E| and show the weight w with the matrix. 

W; = W([ei, e;]) (11) 
And then normalize it to form: 


we W 
~ max(wij EW) (12) 


The network model GN=(Vy, Ey) is initially a regular network that connects the set of edges of Ey 
to each node of the network; although however, this regular network does not offer any good useful features 
for the program. Hence, this regular network must be transformed into a complex network with appropriate 
characteristics. Here we carry this transition using the method applied by Tang ef al. [15] and define the 
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threshold t. This transition is performed to create a new edge set for the network by choosing a subset of 
performed the E*y of Ey which eliminates the need for the edges of the En to have a lower weight than the 
threshold of t. This can be defined as: 


wy =0 ifwj 2t 


Giy = 6(Gy,t =| ; 
CN ( N ) wj=1 if wij <t (13) 


where t varies between tọ and tp and the first and final thresholds (tf, tg) are user-defined. 

Given a graph G=(V, E), we first model the graph structure using network Gy and then extract 
measurements from networks Gy = 6(Gy,t) by varying the threshold t. Afterwards, our CNCRG is 
performed as the concatenation of these measurements. Three distinct feature vectors are proposed, i.e., 
degree descriptors, joint degree descriptors, and clustering-distance descriptors. 

Using formulas described in (1)-(3), @, denote degree descriptors as follows: 


Pt = (Kmax(t), Kavg (t), Ka) (14) 
Using formulas described in (3) and (4), y; denote joint degree descriptors as follows: 

Ve = (Hjalt), Eja(t)) (15) 
Using formulas described in (7) and (8), §& denote clustering-distance as follows: 

& = (C(t), Ga(t)) (16) 
The final feature vector for Giy can be computed as: 


fe = [Pe Ye & (17) 


with the feature vector for Gy, our CNCRG can be computed as the concatenation of f, at different stages of 
the evolution of the network, according to the following: 


CNCRG = fo, fov = fej (18) 


This approach is used to convert the regular network to the complex network, as illustrated in Figure 2. 


0 100 200 300 


Figure 2. The image of a shape and its complex network 


In [14], they showed the complex network Ge is as suitable for the small-world model. This 
complex network model has two main characteristics: i) high clustering coefficient and ii) small-world 
property that can be based on the discussed issues. The clustering coefficients listed in (7) can be used to 
evaluate the characteristics of high clustering coefficients. The small-world property is quantified based on 
whether an average shortest distance is present on the network (8). 
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In this paper, the four algorithms introduced in the previous section for identifying the edge were 
used to obtain the best pattern for identifying the comers and key points, and its output was examined for two 
standard images (camera man and house) in MATLAB software. The results of the implementation of the 
algorithms and their comparison together indicated that the Harris algorithm has higher accuracy and speed, 
and thus, it is used in the rest of the stages. It shows in Table 1. 


Table 1. The comparison of the number of nodes and edges with two standard images 


Number of Edges Number of Nodes Method Seq. 
House Cameraman House Cameraman 
5050 16836 101 184 Harris 1 
19900 27150 200 258 Shi 2 
63903 390286 358 883 SIFT 3 
4186 16110 92 179 SURF 4 


There are several ways to identify neighbors in order to network and connect with them. For 
example, Leandro et al. [25] use the tourist walk method. In this method, first each node is randomly selected 
and then the nearest node is selected based on the Euclidean distance between them. Then the edge between 
these two nodes is drawn. Now for the new node, the nearest node other than the previous node is selected. 
This method continues until they reach one of the pre-selected points. Eventually all of these networks will 
be merged and a new network will be created. Tang et al. [19] used the k-means algorithm to identify 
neighboring nodes. 

This paper used the k-nearest neighbors (K-NN) algorithm to identify neighbors and uses triangles 
to construct triangles; in other words, after identifying the nodes adjacent to one node, the edge was formed 
between the main noun and the neighbors. Finally, after performing this process for all nodes, the parameters 
related to the shortest path, clustering coefficient, and other required parameters were extracted. Figure 3 
shows the steps of the proposed approach. 


Picture 1, Picture2 ,.... PictureN === Edge identification with the Harris algorithm 


4 


Create Complex Network m Identify triangles based on adjacent nodes 


U 


Calculate the properties of the generated graph ==> Classification by neural network 


Figure 3. Classification of images with the proposed hybrid approach 


— Step 1-Edge detection and extraction of key points of the image 

— Step 2-Segmentation of points using Euclidean distance to form triangular areas 

— Step 3-After identifying nodes and their relationship, the network is formed using specified nodes and the 
relationships between them. 

— Step 4-complex network properties extraction (CNCRG) for image description 

— Step 5-The information obtained from the previous step is given for neural network training. For this 
purpose, 70% of information is used for network training and 15% for verification and the remaining 15% 
for testing. 


4. RESULTS 
4.1. Performance evaluation of the proposed method 

System testing is conducted on various sample images available in MATLAB software, for 
example, the corresponding results for two images with their outputs to form a complex network are shown 
in Figure 4 and Figure 5. Assume that our image is a matrix m*n, where m and n indicate the number of rows 
and columns, respectively. For assigning a number to pixel A, which is located in row i and column j, where 
1<i<m and 1<j<n, the number is obtained by the following relation to the pixel: 
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For better comparison, the edges of the image are identified according to Harris and SURF algorithms, which 
part of the results is shown in Tables 2 and 3. 


0 100 200 300 0 100 200 300 


Figure 4. Result of cameraman image processing 


0 100 200 300 0 100 200 300 


Figure 5. Result of home image processing 
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Table 2. Part of sample information with Harris algorithm 
Picture Node No. Triangle No. Clustering Distance Avg. Weight 


1 40 15 3.9 59.86 96 
2 38 15 3.7 58.2 87 
3 43 16 3.97 60.07 97 
4 41 15 3.83 59.98 89 


Table 3. Sample part information with SURF algorithm 
Picture Node No. Triangle No. Clustering Distance Avg. _ Weight 


1 69 22 3.9 72.71 97 
2 62 20 3.7 70.32 93 
3 78 27 3.97 75.43 98 
4 65 23 3.83 71.3 94 


Then information of x and y pixels and their numbers were stored in a table. Then, Euclidean 
distance between them was calculated, and the output was stored in a matrix containing the number of pixels, 
and the distance between them. Then, based on the closest neighboring algorithm, the nodes were identified, 
and practically triangular triangles were formed. Then, the network properties were obtained from the 
formation of graphs based on triangles. 

The above information normalized and output data for each image, including the number of 
identified nodes and triangles, the clustering coefficient, the average shortest distance, and finally, the 
normalized weight were stored in a table. This information was then transferred to the neural network for 
classification. For the initial test, the synthetic sample was first used for testing. An artificial dataset is a set 
of images of a home model using a set of images. Adjacent houses have been achieved with a 2-degree 
rotation. The total number of images was 180 images, with the grid diagram shown in Figure 6. This suggests 
a similar node number but different edge structures of the graphs . 


220 -—— - — 220 -—— -—— ++ 1220 >-———___—_-—— = 
e oa — 
200 - 200 | {200} i 
180 poors 180 | ponosne 1180 | —— 
160 a a 160 m = 1160 | ——— o 
140 v1 140 ci 1140 g E 
120 mul 120 LL 1120} {TL 
100 — n 100 — pof = — 
80 á à à å á á 80 r r à 80 a a a a a 
40 60 80 100 120 140 160 180 40 60 80 100 120 140 160 180 49 6ọ 80 100 120 140 160 180 


100 


40 60 80 100 120 140 160 180 40 60 80 100 120 140 160 180 © 40 60 80 100 120 140 160 180 


Figure 6. Sequence of the model house with feature points in the top row and the associated Delaunay 
diagram with the Harris algorithm 


As for other complex network parameters, the situation was the same, and the case was 
approximately equal to the original image values. The reason for the difference was the rotation of the 
images. To test the system, the coil image set consisted of two classes, and each class consisted of 72 images, 
including duck and box images. The images of these two sets had a 5-degree rotation angle around the 
Y-axis. As previously stated, the table for each image was provided and stored in a table for neural network 
training. 

The results of the duck image with the dual algorithm discussed in Figures 7 and 8 are shown for the 
box image in Figures 9 and 10. The results of the work on these images illustrated the better performance of 
the Harris algorithm. Besides, the similarity for the box-related class was even better due to the box shape 
structure and the clearness of the corners in each rotation, while the information for the duck image was 
distinct, and the number of edges varied in each segment. 
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Figure 7. Results for duck image with harris 
algorithm 
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Figure 8. Results for duck image with SURF 
algorithm 
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Figure 9. Results for box image with Harris 
algorithm 


5. CONCLUSION 
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Figure 10. Results for box image with SURF 
algorithm 


The present study introduced a hybrid method based on complex and neural networks for image 
classification. For this purpose, the key points of the images were first identified. By identifying key points 
using different images of an object at different angles and using the edge-matching algorithm, it identified the 
edges and the formation of triangles. The image is represented as a graph, and by converting the graph into a 
complex grid, the image looks like a complex grid. By converting the image into a complex grid, the 
parameters associated with it can be calculated. Finally, the calculated information was fed to the neural 
network, and it was determined that these parameters are suitable properties for comparing and classifying 
images. It is recommended to use more appropriate approaches and methods for classification by improving 


the neural network learning process or similar methods. 
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